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Abstract. - We study tfie interaction-induced localization - the so-called self-trapping - of a 
neutral impurity atom immersed in a homogeneous Bose-Einstein condensate (BEC). Based on a 
Hartree description of the BEC we show that - unlike repulsive impurities - attractive impurities 
have a singular ground state in 3d and shrink to a point-like state in 2d as the coupling approaches 
a critical value (3* . Moreover, we find that the density of the BEC increases markedly in the vicinity 
of attractive impurities in Id and 2d, which strongly enhances inelastic collisions between atoms 
in the BEC. These collisions result in a loss of BEC atoms and possibly of the localized impurity 
itself. 
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Impurities immersed in liquid helium have proven to be 
a valuable tool for probing the structure and dynamics of 
a Bose-condensed fluid [1]. An example is the use of im- 
purities for the direct visualization of quantized vortices 
in superfluid ^He [2]. Recently, the experimental realiza- 
tion of impurities in a BEC [3, 4] and the possibility to 
produce quantum degenerate atomic mixtures [5, 6] have 
generated renewed interest in the physics of impurities. 
In particular, it was pointed out that single atoms can 
get trapped in the locahzed distortion of the BEC that 
is induced by the impurity-BEC interaction [7-11]. More 
precisely, the impurity becomes self-trapped if its interac- 
tion with the BEC is sufficiently strong to compensate for 
the high kinetic energy of a localized state, an effect akin 
to the self-trapping of impurities in liquid helium [1]. 

However, there are two fundamental differences between 
liquid helium and a BEC. First, for typical experimental 
parameters the healing length of a BEC is three orders of 
magnitude larger than for liquid helium [12], and thus the 
so-called 'bubble' model [13] cannot be applied. A second 
important difference is that the impurity-BEC interactions 



are tunable by an external magnetic field in the vicinity of 
Feshbach resonances [14, 15]. In view of recent experimen- 
tal progress [16] it thus seems possible to investigate the 
self-trapping problem in the same physical system over a 
wide range of interaction strengths - for both attractive 
and repulsive impurities. 

In the present paper we study the self-trapping prop- 
erties of impurities for strong attractive and repulsive 
impurity-BEC coupling within the framework of a Hartree 
description of the BEC. This necessitates the use of an es- 
sentially non-perturbative approach, which is in contrast 
to previous analytical studies [7-9] where the effect of the 
impurities on the BEC was treated as a small perturba- 
tion. We first consider the effect of a highly localized S- 
impurity on the BEC in one dimension. This approach 
indicates that the density of the BEC is substantially in- 
creased in the vicinity of an attractive impurity, which 
enhances inelastic collisions and may result in the loss of 
the impurity atom. In addition, we bring forward a scaling 
argument to show that attractive impurity-BEC interac- 
tions can lead to a point-like ground state of the impurity 
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in 2d and 3d. Specifically, the ground state is singular for 
arbitrarily small attractive coupling in 3d, which seem- 
ingly contradicts the known perturbative results. Finally, 
in order to extend our analytical findings we present nu- 
merical results stemming from the exact solutions for the 
ground state of the system, where unlike in previous nu- 
merical studies [10,11] neither the impurity nor the BEC 
are subject to an external trapping potential. 

Our analysis of the self-trapping problem is based on 
the model suggested by Gross [1,17] describing the impu- 
rity wave-function t) interacting with the condensate 
wave- function ■!/)(x, in the Hartree approximation. The 
model is given by the coupled equations 

= --^VV + «^lxlV + <7l^lV, (1) 
Zmi, 

ihdtx = -^y"-x + ^^m\, (2) 

2ma 

where the impurity-BEC interaction and the repulsive in- 
teraction among the BEC atoms have been approximated 
by the contact potentials kS{x — x') and (?(5(x — x'), re- 
spectively. Here, is the mass of the impurity, mf, is the 
mass of a boson in the BEC and the coupling constants 
K and g > depend on the respective s-wave scattering 
lengths [18,19]. The functions ip{x,t) and xi^^t) are nor- 
malized as 

J dx[V'(x)[2 = iV and J dx[x(x)[2 = 1 , (3) 

where N is number of atoms in the BEC. We tacitly as- 
sume that the properties of the BEC, 6.17. the density 
far away from the impurity, are fixed. The interaction 
strength k on the other hand is considered to be an ad- 
justable parameter, since it may be easily changed in ex- 
periments. 

As a starting point we briefly review and systematically 
extend the perturbative results on the self-trapping prob- 
lem [7~9] based on a variational approach. This method, 
however, yields results that arc of second order in k, and 
hence is not adequate to capture the differences between 
attractive and repulsive impurities. In the main part of 
the paper, we present our non-perturbative and numerical 
results in order to reveal these fundamental differences. 
Finally, we conclude with a discussion of the physical im- 
plications of our findings. 

Perturbative results. — To start we reformulate 
Eqs. (1) and (2) in terms of dimensionless quantities for 
a homogeneous BEC with density uq. The two rele- 
vant length scales in the model are the healing length 
^ = fi./ y^gnomb and the average separation of the con- 
densate atoms s = rtg where d is the dimension 



of the system. Introducing the dimensionless quantities 
X = x/^, tp = ip s'^/'^ and x = X S.'^^^ one obtains the time- 
independent equations 

4, = -iv2vi + /37'lxl'^+lV^I'V^, (4) 
£X = -^y^X + m^X- (5) 

Here, a = mb/ma is the mass ratio, p = n/g \s the rela- 
tive coupling strength, 7 = s/^ and e is the energy of the 
impurity. The energy of the condensate ii^bec, the inter- 
action energy Sint and the kinetic energy of the impurity 
i^kin are given by 

i?bec - 7-'/dx(i[VV'|^-|^|^ + i|^|^), (6) 
Eint = P j dx[x|2|VI', (7) 

i^kin = clx[Vx|^ , (8) 

where all energies are measured in units of gn^ and we 
dropped the tilde for notational convenience, as in the re- 
mainder of the paper. Equations (4) and (5) have a trivial 
solution with = 1 and e — f3, which corresponds to a 
delocalized impurity. 

Thomas-Fermi approximation. If the density of the 
BEC changes smoothly we can apply the Thomas-Fermi 
approximation, i.e. neglect the term V^V in Eq. (4). 
In this regime we immediately find that ['0(x)[^ = 1 — 
/37'^[x(x)[^, and thus obeys the self-focusing nonlin- 
ear Schrodinger equation 

e'x^-^V^x-ClxI'x, (9) 

with the self-trapping parameter ( = and e' = 

(e — P)/a. In one dimension Eq. (9) admits the exact 
solution [7,20] xxi^) = (2A)~i/2gg(,i^(-2,/_)^)^ ^^^^i the lo- 
calization length A = 2/C and the energy e' = — C^/8. 
In addition, numerical solutions of Eq. (9) show that the 
wave-function of the form xa(x) = A^Asech([x[/A), with 
Nx the normalization constant, is also an accurate ap- 
proximation to the exact solution in two and three dimen- 
sions [20,21]. However, these solutions can self-focus and 
become singular in finite time [20]. 

A comparison of the individual terms in i?bcc + Ei^t 
shows that the Thomas-Fermi approximation is valid for 
C "C 1, which is a rather stringent condition on the system 
parameters. However, the merit of Eq. (9) lies in the fact 
that it yields an almost exact wave-function for weakly 
localized impurities with iioc 3> 1, where ^loc is the local- 
ization length of the impurity. In addition, we see from the 
solution of Eq. (9) that self-trapping occurs for arbitrarily 
small coupling in Id for both types of impurities. 
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Fig. 1: The localization length a (solid lines) as a function of 
the self-trapping parameter C, obtained by using a Gaussian 
variational wave-function. Self-trapping occurs for arbitrarily 
small C, in Id, for C, > 2-r in 2d and for C, > 31.7 in 3d (vertical 
dashed lines) . From the two self-trapping solutions for a given 
C, in 3d, only the more localized state is stable. 



Weak-coupling approximation. An alternative ap- 
proach is to linearize the equation for 'ip{'sc) by consider- 
ing small deformations (5'0(x) = V'(x) — 1 of the conden- 
sate [7-9]. This corresponds to a consistent expansion of 
Eqs. (4) and (5) in /3, which results in the linear equations 
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(10) 
(11) 



where we assumed that (5-0 is real for simplicity. It can be 
seen from Eq. (10) that the linearization of Eq. (4) is valid 
in the regime |/3|7''/^foc < 1- The solution of Eq. (10) is 
given in terms of the Green's function G'(x) satisfying the 
Helmholtz equation. One finds by inserting the solution of 
Eq. (10) into Eq. (11) that x obeys the nonlocal nonlinear 
Schrodinger equation 



iv2-2C j dx'G(x-x')|x(x')P 



(12) 



Equation (12) minimizes the energy functional i^[x] = 
E'kin + -Edef , with 



EAoi = dxdx'|x(x)|'G(x - x')|x(x' 



(13) 



which has been discussed in the context of liquid helium 
by Lee and Gunn [7]. We see from Eqs. (8) and (13) that 
the self-trapping parameter C, is proportional to E'def/^'kin, 
i.e. the ratio between the potential energy gained by de- 
forming the BEG and the kinetic energy of the impurity. 

In order to estimate the critical parameters for which 
self-trapping occurs we insert the harmonic oscillator 



X^(x) (ttct^) 



-d/4 



ncxp(-x-^V2a^), (14) 



with a the harmonic oscillator length, as a variational 
wave- function into the functional F[x\. The resulting en- 
ergy /(ct) is given by 



/(a) = ^-/3VM<^): 



with the functions 



/iid(cr) = - exp(2(7^)erfc(\/2cr) , 

h2d{(T) = --^cxp(2a2)Ei(-2(T2) 
zvr 



/l3d(o-) 



27r(T 



exp(2cr2)erfc(V2cr) 



(15) 

(16) 
(17) 

, (18) 



where erfc(a;) is the complementary error function and 
Ei(a;) the exponential integral. The impurity localizes if 
/(cr) has a minimum for a finite value of cr, which depends 
only on the self-trapping parameter Q. 

The explicit expression for /(cr), which appears to be 
a novel result, allows us to determine the self-trapping 
threshold, i.e. the critical coupling strength above which 
self-trapping occurs, and the size of the self-trapped state 
even for highly localized self-trapping solutions as long as 
|/3|7^/o''' ^ 1. Figure (1) shows the localization length 
cr as a function of the self-trapping parameter C in Id, 
2d and 3d. Specifically, in Id we find by asymptotically 
expanding /(cr) in the limit cr ^ 1 that there exists a self- 
trapping solution for arbitrarily small C with cr = jC^. 
The same expansion yields a critical value Ccrit — 27r above 
which self-trapping occurs in 2d. In the 3d case, numer- 
ical minimization of /(ct) shows that the critical value is 
Ccrit — 31.7, where the corresponding state is highly lo- 
calized with a ~ 0.87. Interestingly, there exist two self- 
trapping solutions for a given Q in 3d, however only the 
more localized state is stable. 

Beyond perturbation theory. We now consider 
the self-trapping problem in the regime of strong impurity- 
BEG interactions, i.e. the case of a localized self-trapping 
state accompanied by a possibly strong deformation of the 
BEG. We first illustrate the effect of higher order terms 
in /? for a Id system and subsequently discuss 2d and 3d 
systems based on a scaling argument. 

The S -approximation in Id. In order to find approxi- 
mate solutions of Eqs. (4) and (5) we assume that the im- 
purity is highly localized due to the strong impurity-BEG 
interaction. More precisely, we consider the regime where 
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Fig. 2: (a) Deformation of the BEC for an attractive (left) and 
repulsive (right) (5-impurity with the same interaction strength 
7|/3| = 1. (b) Deformation energy -Edcf for a (5-impurity as a 
function of the interaction strength /3 (solid line) compared to 
the weak-coupling approximation (dashed line) for 7 — 0.5. At- 
tractive (repulsive) impurities have lower (higher) energy than 
predicted by the weak-coupling approximation. 



the effect of the impurity on the BEC is essentially the 
same as for a (5-impurity with |x(cc)P = d{x). This approx- 
imation is valid provided that £ioc <C 1, i.e. the localization 
length is much smaller than the healing length of the BEC, 
which is achieved for sufficiently large |/3|. Our approach 
focuses only on the deformation of the BEC caused by the 
impurity. Nevertheless, the (5-approximation allows us to 
illustrate the differences between attractive and repulsive 
impurities, which are also present in higher dimensions. 

We first calculate the condensate wave- function tlj{x) in 
presence of a (5-impurity at position xq, which is deter- 
mined by the equation 



1 



V(.t)=0, (19) 



with the boundary conditions '4'{x) = 1 and dxi^{x) = in 
the limit x — > ±00. Using the known solutions of Eq. (19) 
for (3 = and taking xo = for simplicity we obtain 
il}{x) = coth(|a;|-|-c) and tpi^) = tanh(|a;|-|-c) for attractive 
and repulsive impurities, respectively, with c a constant. 
Given that the derivative of il^{x) has a discontinuity at xq 
it follows that in both cases 1 — [^'(xo)]^ = l3^tp{xo) and 
thus 



If ^ 

V[xo) = — ^ 



1 



(20) 



which fixes the constant c. We note that ^'(2^0) G (0, 1] for 
/3 > and ipixo) G (1, 00) for /3 < 0. Thus the density of 
the condensate n{xo) = [V'(xo)]^ at the position of the im- 
purity xq is not bounded for attractive interactions. The 
enhanced deformation of the BEC due to an attractive 
(5-impurity is already visible for a moderate interaction 
strength as illustrated in Fig. (2)a. 

The solution of Eq. (19) allows us to determine the de- 
formation energy E'dcf to all orders of /?. Using Eqs. (6) 



and (7) we find that for both types of impurities 

3/2 



(21) 



We note that the non-perturbative result for Edof is con- 
sistent with the weak-coupling approximation since in the 
limit (T ^ we find from Eq. (16) that (3'^jhia{a) 
/3^7/2. Figure (2)b shows the deformation energy _Edcf 
as a function of the interaction strength /3 in comparison 
with the weak-coupling result. It can be seen that attrac- 
tive impurities have lower energy than predicted by the 
weak-coupling approximation, whereas repulsive impuri- 
ties have higher energy. 

Impurity collapse. The results in the previous section 
indicate that the ground state of the impurity and the 
BEC strongly depend on the sign of /?. We now show that 
the sign of /3 has a crucial effect on the ground state in 2d 
and 3d by using a scaling argument based on a variational 
wave-function for 'ip{x.) and x(x). 

Let us consider the total energy of the system _Etot = 
E'bcc + -Eint + -Ekin , which is clearly bounded from below for 
(3 > but not necessarily for (3 < 0. We now analyze the 
scaling of the individual terms in i^tot depending on the 
localization length of the impurity and the deformation of 
the BEC. To this end we insert the Gaussian trial function 
Xct(x) for the impurity and V'(t(x) = 1 + (5^/v(x) for the 
condensate into -Etot; where 



(5?/>o.(x) 



rS/2 



Y[ eyLpi-x^ba^) , 



(22) 



with a, b and S positive constants. In what follows we 
are particularly interested in the limit cr — > and only 
consider finite deformations of the BEC. Consequently, we 
add the constraint 6 < d to assure that /dx|(5'0(x)p 
(t'*^* is finite. We note that Xct(x) and V'o-Cx) have the 
correct asymptotic behavior required for them to be valid 
physical states of the system, however, they may not be a 
good approximation for the ground state. We obtain that 
the individual energy contributions scale as Ei^t ~ cr^^ , 
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kin 



and 
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J = l 



d-]S/2 



(23) 



where the Cj > depend on the system parameters. 

Comparing the terms above one finds that for a three- 
dimensional system i?tot is not bounded from below in 
the limit ct — > for a fixed 6, i.e. the energy of the 
system becomes arbitrarily low if the impurity collapses 
to a point. Explicitly, in the 3d case we have -Ebcc ~ 
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Fig. 3: The localization length a (solid line) and the density of 
the condensate n{xo) (long-dashed line) at the impurity posi- 
tion xo as a function of the interaction strength f) in Id. The 
weak-coupling approximation (dotted line) accurately repro- 
duces cr in the regime |/3|7/(7 ^ 1. The localization length 
saturates to a finite value cr ~ 60 as /3 — ^ because of the 
finite system size. The deformation of the BEC is qualitatively 
described by (5-approximation (short-dashed line). The inset 
shows the smooth transition from a sech-type {Rm = 1) to a 
Gaussian [Rat = — 1) self-trapping state for repulsive impuri- 
ties. 



Fig. 4; The localization length a (solid line) and the density 
of the condensate n{xo) (dashed line) at the impurity position 
Xo as a function of the interaction strength (3 in 2d. The weak- 
coupling approximation (dotted lines) accurately predicts the 
self-trapping threshold at /3crit ~ V^Stt (vertical dashed lines). 
The localization length saturates to a finite value cr ~ 20 for 
\P\ < /3crit because of the finite system size. Attractive impu- 
rities shrink to a point-like state and the density of the BEC 
n{xo) diverges as /3 approaches a critical value f3* ~ —10 in 
agreement with our scaling argument. The behavior of Rat 
(inset) is similar to the Id case. 



cq(t^~^ + cia^^^/^ -I- • • • -I- C4cr'^~^'' and hence a compari- 
son of the exponents yields that ji^intl > -E-kin + E^oc for 
2 < (5 < 3 in the limit a ^ 0. Moreover, in the same limit 
J dx\6ipc{x)\'^ 0, i.e. the deformation of the BEC is 
vanishingly small despite the collapse of the impurity. We 
note that the above argument holds for arbitrarily small /3, 
and thus the weak-coupling approximation in 3d is valid 
for /3 > only since for (3 < there is always a state of 
lower energy with a ^ 0. 

In the 2d case it is not possible to find a value for S 
such that the total energy i?tot is dominated by Eint in 
the limit cr ^ 0. Nevertheless, for (5 = 2 we obtain i^int ~ 
-Ekin ~ £'bec ^ cr^^, which implies that the collapse of the 
impurity depends on the particular system parameters, 
e.g. the couphng strength (3. In contrast to the 3d case 
this collapse is accompanied by a sizeable deformation of 
the BEC since / dx|(5'(/'cr(x)P is finite in the limit ct ^ 
for 6 = 2. Finally, for a one-dimensional system similar 
considerations lead to result that |i?int| < ^^kin + -Ebcc in 
the limit cr ^ for any valid S, and hence our variational 
approach docs not predict a collapse of the impurity. 

Numerical results. — We have determined the 
ground state of the impurity and the BEC numerically 
by using the normalized gradient flow method [22, 23] in 
order to extend the analytical results. To make Eqs. (4) 
and (5) amenable to computational modelling we assumed 
that both the impurity and the condensate were enclosed 



in a sphere of zero potential with infinitely high poten- 
tial walls at its radius R. The system parameters for the 
numerical calculations were set to a = 1.0 and 7 = 0.5, 
whereas (3 was varied over a large range. 

A quantitative measure for the localization length ^loc 
was found by fitting the fimctions Xct(x) and Xa(x) to 
the exact impurity wave function yielding the values a 
and A. In addition, we compared the corresponding fit- 
ting errors So- and S\ and evaluated the quantity i?fit = 
{Sa'~S\)/{Sa + Sx). Accordingly, i?fit = —1 indicates that 
the self-trapped state is accurately described by a Gaus- 
sian X(t(x), whereas for i?fit = 1 the impurity is rather in 
a sech-type state xa(x). One would expect that i?fit « — 1 
for a highly localized impurity with £ioc <C 1 since the 
impurity sees only the harmonic part of the effective po- 
tential /?|7/'(x)p. On the other hand, Rfn « —1 for a 
weakly localized impurity with £ioc ^ 1 according to the 
Thomas-Fermi approximation. 

As can be seen in Fig. (3), in the Id case there is 
no self-trapping threshold in agreement with the weak- 
coupling result, which accurately reproduces the localiza- 
tion length cr in the regime |/3|7/cr <c 1. However, we see 
that for large attractive impurities are more localized 
than the repulsive ones. Moreover, repulsive impurities 
undergo a smooth transition from a sech-type to a Gaus- 
sian self-trapping state for increasing (3 as indicated by 
i?fit. The deformation of the condensate n(xa) at the po- 
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Fig. 5: The localization length a (solid line) and the density 
of the condensate n{xo) (dashed line) at the impurity position 
Xo as a function of the interaction strength j3 in 3d. There 
is no well-defined ground state for /3 < in 3d in agreement 
with our scaling argument. The weak-coupling approximation 
(dotted line) predicts a self-trapping threshold at /3crit — 15.9 
(vertical dashed line) , which is lower than the numerical result 
/3crit — 18.8. The localization length saturates to a finite value 
CT ~ 10 for /3 < /3crit because of the finite system size. The 
convergence of the impurity wave-function to a Gaussian self- 
trapping state is slower than in Id and 2d as indicated by Rm 
(inset). 



sition of the impurity xq is qualitatively described by the 
(5-approximation with a strong increase in the density of 
the BEC for attractive impurities. 

In two dimensions, the impurity localizes for > /3crit, 
where Pent is correctly predicted by the weak-coupling re- 
sult, as shown in Fig. (4). This is not accidental since the 
localization length cr diverges close to Pcrit in 2d, as can be 
seen in Fig. (1), and thus the criterion |/3|7^/fT^ ^ 1 for 
the validity of the linearization of Eqs. (4) and (5) is always 
met near /3crit- In contrast, for large \(3\ the localization 
length (T strongly deviates from the weak-coupling result. 
In particular, attractive impurities shrink to a point-like 
state and the density of the condensate ri(xo) diverges as 
/3 approaches a critical value (3*, which is in agreement 
with the scaling argument in the previous section. The 
behavior of Rat is similar to the Id case. 

For a 3d system we are able to obtain numerical results 
for /3 > only, which in view of our scaling argument was 
to be expected. The localization length a and the density 
of the condensate n(xo) for repulsive impurities are shown 
in Fig. (5). We see that the weak-coupling approximation 
underestimates the exact value /3crit since in 3d we have 
that CT ~ 1 even close to /3crit and thus P^^/a^ is not 
necessarily small. As shown in the inset of Fig. (5) the 
transition from a sech-type to a Gaussian self-trapping 
state is slower than in Id and 2d. 



Discussion and conclusion. In our analytical and 
numerical investigation of the self-trapping problem we 
have shown, in particular, that an attractive impurity- 
BEC interaction leads to a strong deformation of the BEC 
in Id and 2d and to a singular ground state of the impurity 
in 2d and 3d. However, both the high density of the BEC 
and the point-like state of the impurity entail additional 
physical effects, which were initially small and not taken 
into account in our model. 

The high density of the BEC near the impurity enhances 
inelastic collisions that lead to two- and three-body losses 
of the condensate atoms [24]. Specifically, the inelastic 
collisions might cause the loss of the impurity atom it- 
self. However, provided that the impurity is not affected 
by inelastic collisions one can account for the loss of the 
condensate atoms by adding damping terms of the form 
-ir2|-0(x)|V(x) and -ir3|?/>(x)|'^-0(x) to Eq. (1), where 
the positive constants r2 and are two- and three-body 
loss rates, respectively. We note that these additional 
terms would require the simulation of the full dynamics 
of the system [25] and may have a nontrivial effect on the 
self-trapping process. 

Moreover, the contact potential approximation is valid 
as long as the various length scales of the system, e.g. the 
localization length ^loc, are large compared to the charac- 
teristic range of the exact interaction potential [18]. Thus 
the interaction between the condensate atoms and the im- 
purity is no longer correctly described by the contact po- 
tential K(5(x — x') in the limit iioc — > where the influence 
of higher partial waves becomes important. The resulting 
repulsion would of course lead to a finite width of the im- 
purity wave- function. We note that with the above proviso 
our results are also valid for ionic impurities with positive 
charge [26,27] in the regime where the ion-boson scattering 
is dominated by the s-wave contribution [28,29]. 

Finally, we conclude with the remark that our model 
is also applicable if the BEC is confined to a harmonic 
trap [10] provided that the harmonic oscillator length of 
the potential is much larger than the localization length 
£\oc of the impurity and the healing length ^ of the BEC. 

^ ^ ^ 
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